Spin-Filtering Multiferroic-Semiconductor Heterojunctions 



O 

o 
c/3 



o 

CZ3 



o 
o 



> 

oo 

m 
p 

ON 

o 
o 



X 




Na Sai, 1 Jaekwang Lee, 1 Craig J. Fennie, 2 and Alexander A. Demkov 1 

1 Department of Physics, The University of Texas at Austin, Texas, 78712 
2 Center for Nanoscale Materials, Argonne National Laboratory, Argonne, IL, 60439 

(Dated: August 24, 2007) 

We report on the structural and electronic properties of the interface between the multiferoic 
oxide YMnC>3 and wide band-gap semiconductor GaN studied with the Hubbard-corrected local 
spin density approximation (LSDA+U) to density- functional theory (DFT). We find that the band 
offsets at the interface between antiferromagnetically ordered YMnC>3 and GaN are different for 
spin-up and spin-down states. This behavior is due to the spin splitting of the valence band induced 
by the interface. The energy barrier depends on the relative orientation of the electric polarization 
with respect to the polarization direction of the GaN substrate suggesting an opportunity to create 
magnetic tunnel junctions in this materials system. 



Multiferroics are materials in which ferroelectricity and 
magnetism coexist in a single phase. Efforts have shifted 
from the question of "coexistence" , to identifying strate- 
gies that will increase the coupling between the two or- 
derings. This increased functionality opens opportunities 
for novel electrically or magnetically controlled devices. 
One route that is promising for practical applications is 
to employ multiferroic thin films and multilayered struc- 
tures, the properties of which can be readily manipulated 
at the nanoscalei 2 - Multiferroic YMnC>3 3,4 is considered 
an attractive candidate for use in transistor devices be- 
cause of its purity from volatiles such as Pb or Bi and its 
moderate dielectric constants. 

Much effort has been directed at synthesizing and char- 
acterizing thin films of hexagonal YMnC>3 as a potential 
gate dielectric for semiconductor devices, e.g.,YMnC>3 on 
Si(lll)£&I£, on wurtzite GaN, or on ZnC&I&II. De- 
spite the remarkable progress in synthesis, the role of 
the interface — strain, chemistry, etc., — on the ferro- 
electric and magnetic properties of YMn03 thin films is 
poorly understood and difficult to separate experimen- 
tally. In this Letter, we apply density-functional theory 
to calculate the electronic structure and band alignment 
at a realistic YMnOs/GaN interface. We demonstrate 
that interfacial spins behave differently from those in the 
bulk. This interface effect leads to a spin splitting in 
the valence bands giving rise to different band offsets for 
spin up and down states. Intriguingly, the difference in 
the band offsets depends on the polarization direction of 
YM11O3 relative to that of the polar GaN substrates, sug- 
gesting that the system could be utilized in spin-filtering 
tunneling junctions. 

Bulk YMn03 and wurtzite GaN both have hexagonal 
symmetry. YMn03 is antiferromagnetic (AFM) and fer- 
roelectric (FE) (space group P63CTO) while GaN is polar 
but not FE. X-ray diffraction of YMn03 thin films de- 
posited on GaN (0001) substrates suggests an in-plane 
rotation of 30° between the unit cell axes of YMnC-3 and 
GaN. For a coherent interface this implies that YMn03 is 
under a 10% in-plane compressive strain^ ' 10 ^ 1 This large 
epitaxial strain has been attributed by Posadas et al. to 
the energy gain from the interfacial bond formation. 11 
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FIG. 1: The structural model of the YMn0 3 -GaN heterojunc- 
tion. The two inequivalent interfaces viewed along (0001) and 
(0001) are defined as A and B, respectively. 



Here, we adopt the experimentally determined interfa- 
cial relation and build heterostructures composed of two 
unit cells of YMnC-3 and GaN each, with a total of 24 
atomic layers (Fig. [IJ. Both YMnO^/GaN interfaces 
are comprised of Ga and apical oxygens, coming from 
YMnOs and GaN, respectively. The interfacial oxygens 
are placed above the fee sites (i.e., the threefold cavity 
sites) of the Ga surface. The periodic boundary condition 
in our calculations creates two inequivalent Ga-0 bonded 
interfaces between the sequence of YM11O3 and GaN, one 
with oxygens above the Ga (0001) face and another with 
oxygens above the Ga (0001) face (see interfaces A and 
B in FigQ] respectively). 

We perform DFT calculations of YMn03 within the 
LSDA+U approximation with U = 6 eV and J = 0.9 
eVi^ 2 - The value of U was extracted from experiment 
(XPS^ 3 -). We obtain a band gap of 1.46 eV (1.47 eV^ 3 - 
and lattice parameters of oymo = 

6.09A (6.127A)±£ and 
c/a = 1.86 (1.86) which agree well with experiment 
(shown in parentheses). The Mn spins were treated in a 
frustrated collinear-AFM approximation as the 120° non- 
collinear spin structure observed in bulk YMn03 is be- 
yond the computational capabilities for the realistic inter- 
face that we consider here. For bulk GaN, we carried out 
LDA calculations and found acaN = 3.15A, c/a = 1.627, 
and u = 0.377. The LDA band gap is 2.1 eV, which is 
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FIG. 2: The profile of the macroscopically averaged electro- 
static potential across the YMn03-GaN supercell Pi. The 
straight lines are the top of the valence bands measured with 
respect to the average of the potential in each material, the 
values are referenced from independent bulk calculations. 



below the true value 3.5 eV. Since the Hubbard correc- 
tion does not provide genuine improvement in GaN, we 
used a Hubbard U only on the Mn atom. 

To isolate the effect of strain, we first relax 
bulk YMnC>3 under in-plane compressive strain of 
VScLGaN I flYMnO ~ 10%. The c-axis lattice parameter 
expands by ~7%, consistent with Ref. |Tl|. The Mn mo- 
ments are slightly reduced ~ 3.6^p compared to the un- 
strained bulk value of 3.77/ip. The distances between Mn 
and the in-plane oxygens Mn— Op shorten to ~ 1.86A, 
while the distances between Mn and the apical oxy- 
gens Mn— Oa lengthen to ~ 1.9A, compared to ~ 2.04A 
and ~ 1.86A, respectively, in the unstrained bulk. De- 
spite the significant shortening of the Mn-Op bonds, no 
changes occur to the topmost valence bands including 
Op Px,Py and Mn d xy and d x 2_ y 2 orbitals as well as the 
unoccupied d z 2 states. The most noticeable change oc- 
curs in the O^ p states and the Op p z states immediately 
below the Op p x<y bands, which downshift by ~ 0.6 eV 
near the T point, presumably due to strain- induced buck- 
ling of Mn05 bipyramids. 

Because wurtzite GaN is polar (a permanent dipole 
moment is produced by each bilayer of Ga + and N — ) we 
construct two supercells; "Pi" with the FE polarization 
in YM11O3 pointing along GaN (000T), and "P 2 " with the 
FE polarization pointing along GaN (0001). We then re- 
lax both supercells fixing the in-plane lattice constant to 
a GaN • In either supercell, the vertical distances between 
the apical oxygens and the Ga atoms at the A and B in- 
terfaces (see Fig(TJ are drastically different, ranging from 
1.1-1.2A at interface A to 0.2-0.6A at interface B. The 
different bonding structures of two interfaces agrees with 
the behavior of oxygen adatoms absorbed on Ga surfaces 
in the low oxygen coverage (growth) condition of GaN^ 

Fig. [5] shows the macroscopically averaged electro- 




FIG. 3: Orbital-resolved DOS projected on atomic layers of 
the YMn03-GaN supercell shown in Fig.Q] In each panel, the 
top and bottom most graph corresponds to interface A and 
B, respectively. The projected DOS of Mn d xy and d yz (not 
shown) are nearly degenerate with d x 2 and d xz respectively 
shown in d). The vertical lines mark the Fermi level. 



static potential the z direction. The polarization field 
in GaN, evident from the slope of the electrostatic po- 
tential, induces charge transfer across the film and in- 
duces screening on the opposite sides of the GaN films. 
As a consequence, the residual field across YMn03, par- 
allel (antiparallel) to the field in GaN in supercell Pi 
(P2), is small. Using the effective ionic charges^ we es- 
timate the polarizations of YMn03 in the supercell to 
be -9.99^C/cm 2 for Pi and 8.43^C/cm 2 for P 2 . These 
numbers are slightly lower than 14.0/LtC/cm 2 calculated 
for the constrained bulk (P ~ 8.8/^C/cm 2 in uncon- 
strained bulk)^, indicating incomplete screening of the 
polarization charges at the interfaces. The polarization 
fields in these systems complicate comparison of the band 
edges directly from the lineup of the average potential as 
shown in Fig. [5] 

To analyze the electronic structure of the supercells 
we calculate the orbital-resolved, layer-projected density 
of states (PDOS). In the Pi supercell, there is clearly a 
band bending of ~ 1 eV in GaN due to polarization, as 
illustrated by the deep Ga id valence states in Fig. |3Jz. 
The band bending is also evident in the topmost valence 
bands N p in Fig. 0> that produces an upward bending 
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FIG. 4: Projected DOS of the YMn0 3 -GaN Pi supercell into 
4 bilayers of GaN (red) and 2 unit-cells of YMnOs (blue). 
The left and right panels are DOS for spin-up and spin-down 
components. The black dots mark the valence band edge and 
the arrows mark the interface band edge. 

at interface A and a downward bending at interface B, 
consistent with the experiments in GaN (0001) filmsjil 
Despite this bending, there is no sign of gap closing as 
the thickness of our GaN films is much below critical^ 
Fig. [3J; and d show the DOS for Mn 3d and planar Op 
2p states. The interior layers look similar to the bulk 
layers, suggesting that the interior region converges to 
the bulk. The hole states at interface A are composed of 
Mn d xy , d} c 2_ y 2, Op p\ y and a small contribution from 
Oa P x y At interface B we observe a down shift in energy 
of the Mn 3dz 2 orbitals (located at ~2eV above Ep in 
the bulk), and a strong overlap with the O p states. We 
encounter similar behavior in supercell Pi except that the 
Ep is ~ 0.2 eV deeper than in PI. The difference can 
be related to charge compensation at the interfaces: in 
Pi , the polarizations of YMn03 and GaN create opposite 
screening charges that partially cancel at the interfaces, 
while in P2 , the polarizations create screening charges of 
the same sign and thus lead to higher accumulation than 
in Pi. 

In both supercells we observe that the spin moments 
of Mn at the YMn03 / GaN interfaces deviate from their 
bulk values. On average the moment of Mn reduces to 
~ 3.3/xp at interface A and increases to ~ at inter- 

face B, producing a small net magnetization of 1.5fip per 
supercell. This can be understood from the DOS analy- 
sis which shows a depletion of charge on the Mn d xy and 
d x 2 at interface A and accumulation of charges on the 
Mn d z 2 orbitals at interface B. It is certainly possible 
that the precise route to the observed ferrimagnetic mo- 
ment depends on the collinear approximation used in this 
study, but the basic physics of spin manipulation at the 
YMnOs/GaN interface should be qualitatively similar if 
noncollinear spins are considered. 

To calculate the band offset at the YMn03 / GaN inter- 



face, we project the spin-resolved DOS of the supercell 
onto unit cells of YMn03 and bilayers of the GaN as 
shown in Fig. 0] and calculated the band offsets for both 
the spin-up and spin-down components. In the case of 
Pi, the valence band offset (VBO) is —0.05 eV for the 
spin-up component and —0.1 eV for the spin-down com- 
ponent at interface A, whereas a VBO is +0.35 eV for 
the spin-up and +0.15 eV for the spin-down at inter- 
face B ("+" and "— " represent upward and downward 
offsets going from GaN to YMn03). The band offset at 
both interfaces shows spin dependence. In particular, the 
spin-up electrons experience a potential higher by 0.2 eV 
than the spin-down electrons across interface B. This 
suggests that electrons of different spin directions expe- 
rience different electrostatic scattering probabilities and 
that spin- dependent barriers can be produced in AFM 
multiferroic/semiconductor heterostructures. This effect 
is associated with spin splitting at the top of the valence 
bands. Recently, BiMnOs tunnel barriers were shown to 
be potential spin filters in magnetic tunnel junctions. 19 
But FM multiferroics such as BiMnOs are rare. Thus our 
observation of spin-dependent barriers in AFM multifer- 
roic/semiconductor junctions might enable the use of a 
larger group of multiferroics with AFM orderings as spin 
filtering devices. 

In the supercell Pi (not shown), we find a VBO of 0.1 
eV for the spin- up and —0.1 eV for the spin-down at in- 
terface A, similar to that of Pi. At interfaces B, the 
VBO is +0.2 eV for spin-up and +0.4 eV for spin-down, 
slightly higher in amplitude than that of Pi. It is note- 
worthy that the offset of the Pi structure is reversed for 
up and down spins from that in Pi at interface B, which 
suggests a possibility of controlling the spin scattering 
by a reversal of the ferroelectric polarization in YMn03. 
A small energy difference of ~ 20 meV is found in bulk 
YMO3 when the spins are switched from the collinear 
to noncollinear in-planc configurations. 2 - This is much 
less than the barrier differences we found between the up 
and down spins. Thus the magnetic anisotropy should 
not influence the spin-dependent effect we observe. 

In conclusion, we report the details of the electronic 
structure of hexagonal YMn03-GaN hetero junctions. 
We consider two inequivalent Ga-0 terminated interfaces 
that can be found in YMn03 films grown on (0001) and 
(0001) oriented GaN substrates and two possible orien- 
tations of the YMnOs polarization with respect to that 
of the GaN substrate. We find different band offsets for 
spin-up and spin-down components, with a larger vari- 
ance at the (0001) interface. The spin-dependent inter- 
face barriers suggest that these heterostructures may be 
applicable in spin filtering tunneling devices. Our results 
are relevant not only to YMnO films but also to other 
multiferroic thin films with coexisting antifcrromagnetic 
and ferroelectric structures. 

This work is supported by the Office of Naval Re- 
search under grant N000 14-06-1-0362 and Texas Ad- 
vanced Computing Center. 
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